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ieee ODUCT ION 


irewrieldwol Signal Processing is expanding rapidly, 
mcdseating an increasing importance. Today, digital signal 
processing is being extensively explored and employed. 

High speed digital computers with necessary peripheral 
equipment such as analog-to-digital and digital-to-analog 
converters are required. Traditional analog systems as radar 
and sonar are going digital. There exists a wealth of other 
digital signal processing applications. 

For some applications, however, the use of digital 
processing may be termed "overkill", even if there is a 
distinct requirement for the inherent advantages of digital 
processing over analog processing. 

A fairly new device technology, the Charge Transport 
Device technology (CTD's) offers new ways of signal processing 
which could be the answer for such cases. 

Charge Transport Devices refer to several families of devices: 
Bucket Brigade Devices (BBD's), Charge Coupled Devices (CCD's), 
and Reticon Devices (i.e. SAD-100 and TAD-12) which are 
described in reference {3}. 

The unique feature of CTD's is that the signal is 
represented by discrete charge packets, which are moved 
through the device at a speed determined by the clock 
frequency. The delay of the device is thus controlled by the 
clock frequency, and may be easily and momentarily changed. 
An excellent review of how the CCD works is given in 
reference {8}. 

CTD's fing wide applications in imaging (both low light, 
normal light and infrared), analog time delay, multiplexing, 
transversal filtering, with fixed and variable weighting, 
recursive filtering, correlators, chirp-z transformers and 
SO on. The general advantages of CTD signal processors 
over conventional digital processing are lower cost, power, 


weight and size. 
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Comb filters are described by their periodic transfer 
eiiearacterasulcs an the frequency domain. They can be classi- 
fied into two general types. The first type is the bandstcp 
Semeanceiicer type shown in Table II1. It has strong attenu- 
ation in a narrow neighbourhood of a series of periodically 
separated frequencies and good transmission in between. 

The second type is the bandpass, or integrator type (Table II). 
It has good transmission in a narrow neighbourhood of a 
series of periodically separated frequencies and strong 
attenuation in between. All possible implementation 
techniques for such comb filters use some kind of delay de- 
VaLCe . 

In the analog case, a quartz delay line has been 
extensively used. It is, however, temperature dependent 
and results have in general not been very satisfactory. 

For the digital case the delay device is the shift register, 
in which the amount of delay is determined by the employed 
clock frequency. 

Today a new family of comb filters is being developed, using 
the new CTD's. In these, analog signals are first sampled 
and then delayed. As a result of the sampling, these devices 
not only have standard analcg properties but also some 
digital properties such as aliasing and stability of delay. 
A CTD consists of N serial delay stages. A signal will then 
be delayed N clock periods before it is fed back to the 
input in the recursive case. The frequency of recursive 
operation, Jas 1s then only Set where iE is the clock 
frequency, or sampling frequency. Therefore, in radar terms, 
on is equal to the Pulse Repetition Frequency (PRF), and 
there are N teeth of the comb within the Nyquist frequency 


range from to [/2. 


Pemmeapemecursave filtering using CCD's have not been 
extensively studied in the industry, since it was believed 


that it can only be done with extensive external electronic 
support such as fast and accurate operational amplitiaers ete. 
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In fact, reference {4} was the only paper on CCD- recursive 
filters presented at the recent conference on CCD-appli- 
cations in San Diego. 

This thesis will study the second order recursive comb filter 
and employ a simulation scheme to study its application 


as a canceller in an MTI-radar. 





I~. SAMPLED ANALOG SECOND ORDER RECURSIVE COMB FILTERS 


A. ie OLY 


memernalee ri lter 


Since the sampled analog recursive filters are 
impopemented by the Same circuit configuration used for digital 
recursive filters, we may apply the existing well~developed 
theory for such filters. However, proper modifications must 
be made to account for the fact that the signal is now sampled 
analog and that there now are N delay stages instead of one. 
It is known that two general design approaches have been 
developed for the IIR (infinite impulse response) type digital 
recursive filter, namely the indirect and direct methods. 

The indirect method starts with the transfer function in the 
LaPlace transform variable s (H(s)), of a continuous 
filter and follows by performing another transform to obtain 
the corresponding transfer function in the discrete time 
variable z. Several different transforms are available 
for this purpose, for example mapping by differential 
transform, impulse invariant (standard z ) transform, 
bilinear z-transform and matched z-transform. 
The direct method is digital design in frequency or time 
domain using some type of computer aided design procedure. 
Here the indirect approach will be used: 
Given a prescribed response, there are several methods by 
which one can develop a transfer function in the variable s. 
Examples are: Butterworth, Chebyshev and elliptical filters 
eeewiinelma this thesis the Butterworth approach is used to 
obtain the transfer function in the s-domain. 
a. For a Low Pass Fiiter 
How the Butterworth filter im the second order 


case the equation is given as: 





Gt a 
aL 





where w, is the 3-db cutoff frequency. 


Converting to s-domain and expanding: 


MCs mi (=s)- = i 


2 N cert N = 2 


tos 
a 
By evaluating the poles in the s-plane and discarding the 


two poles in the right hand plane, one gets: 


HiCso a= I. 
5 a IN Sap 2 
x x 
This is the desired transfer function in the s-domain. 
b. For High Pass Filter 
A similar derivation for the high pass case 


gives: 


a. Transformation to Z-domain 


a. Low Pass Filter 
By indirect design method, using the bilinear 


z-transform: 


f 8 Ue 
TCz+1) 
one gets the following transfer function in 2: 
2 
ae 
IU ¥, 
2 
Oa Ivar Sh ae 
253) * FF zr oy * 


By employing a little algebra one will end up with this 


equation: 
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where: 


K ( a_ + 1 z + ay z ) 
-~1 -2 
1 + bo We + bo Zz 
T? W 2 
X 
u+ 2/2 w T + Ww? ~? 
a 
2 
AL 
2 WwW Q 7? 8 
x - 
i) oy ane: Ww’ Tt? 
4 + a tT? - 2vDwT 
x x 
Dee 


ee ¥2w_T + WwW ut 


ey 
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b. High Pass Filter 


Starting from the s-domain equation: 


Z 
S 


Eis) = = yp te 
Ba 
Ss + ¥2 wo, S +o, 


and using the bilinear transform again, one end up with: 


K ¢ SL ar ay zt + ay Bea 
B64 ees 
-i -2 
3 ee by Z. + be Z, 
where: 
K = = 
i DW aie ee yee 
x se 
ay = an ae 
ay = -2 
2 wT - 8 
by = 


Lian ey wT er 


+ 2 VS oy Ox ey ak 
x x 


ies noted that WW, and T always occur together in the final 
expressions and that there are only two differences between 
the high pass and low pass cases: The size of K, and the slgn 
of ay: But K is merely a scaling factor, and has no influence 
upon the shape of the frequency response. 


This means that when the coefficients are kept constant, the 
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Region II 


Region I 
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Sweecrtearequency wail change by varying T. 


Mamole - and Fagure 2.1 shows how by and be varies with wT. 


Table II is a summary of how the frequency response behave 
as wT is varied, giving the ranges where the filter acts as 
integrator or canceller. 

Beuee region I (Figure 2.1) is the area of interest in this 


thesis, the experimental study was carried out here. 
ee mect | Pesien Method in Z-Domain 


a. By Trial and Error Method 
By carrying out the design in the z-domain, 


feemstormations from H(s) to H(z) can be avoided: 





PuGU beeen 2. ee Z-Domaim Unit Circle 


Figure 2.2 shows the z-plane, where the entire left half of 
the s-plane has been transformed into the inside of the unit 
circle. By moving the two poles and two zeros (second order) 
around and using the geometric method {1} to calculate the 
amplitude- frequency response, it is possible to find a 
combination of the singularities that gives a response quite 
close to the- desired response. Appendix A shows the result 


Sessuch a study. 


b. By Computer Approach 
Several computer programs {1} have been 


developed that will optimize the search of coefficients, 
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VEE eT 


Relationship between wT and Filtertype for Highpass 


and lowpass filters using the bilinear z-transform method. 


Filter Type Ce avilietete eletshare. 
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(Region I, Integrator) 
LOW Pass 
i 
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(Region II, Canceller) 
i | 
ae 


Ee = sampling frequency N= number of delay stages 
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based on different approaches in reducing the deviation 
between presented response and filter response. 


This approach was not employed here. 


B. EXPERIMENT 


le system Configuration 








FIGURE 2.3. Second Order Recursive Filter 


Figure 2.3 shows the block diagram of the canonical form of 
a second order recursive filter. Serial analog delay (SAD) 
from Reticcen Co were used as delay lines. 20 K potentio- 
meters were used to implement the coefficients. In the cases 
where coefficients greater than one were needed for the "a"- 
See*ttctents, the numerator in the transfer function H(z) was 
scaled down so that the "a"-coefficients were normalized. 
several methods to implement the coefficients using the 
potentiometers are discussed in references {2} and {3}. 

In this work Figure 2.4 is presented, giving the calibration 
of potentiometer setting for desired coefficient value. This 


curve was obtained by simply measuring input and output 


Bs 
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across the potentiometer. 
2. Results 


iiesctock frequency for the SAD-100 (Reticon Co) 
device was 72 KHZ. Since the effective clock frequency was 
meurs times 72 KiZ, and there was 96 stages within the device, 
the corresponding nulls of a high pass comb filter should 


mecur at multiples of: 


(2 
— oa CF 3 KHZ 


A computer program (APPENDIX B) was used to calculate and 
plot the theoretical frequency response curve,and also to 
superimpose the actual data points on the same plot. 

The "a"- coefficients were kept fixed, and several combina- 
rons. Of b, and b, were used. They were selected from 
Peewee 2,1, and Dy was negative in all cases studied in this 
thesis. 

By using the values of by and b, corresponding to wT, the 
resulting response should be flat-topped in every case. 


The implemented coefficients were: 
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ime coefficients b, and by Vege c lp leked am) hepaon | 
Meatcure 2.1). 

Fagures 2.9 through 2.18 show the series of results. 

There seemed to be the same problem with the second order 
filter as was experienced with the first order filter{2}, 
namely a skewness at higher frequencies. Also the results 
Show an undesired dip in the passbands, indicating that the 
coefficients are not correctly implemented. 

Suspecting that the experimental implementation inherently 
makes the system frequency dependent, several checks were 
made for its investigation. The gain of individual delay lines 
was checked for frequency dependence, using both the oscillo- 
scope and a narrowband voltmeter. 

The individual potentiometers input-output characteristics 
were measured when the filter loop was closed. 

It was found that the gains of the delay lines were slightly 
frequency dependent as measured on the oscillosccpe, and 
severely and inconsistently frequency dependent when 
measured on the narrowband voltmeter, indicating loading 
effects in the latter case; (the input impedance of the 
voltmeter was 100 Kohm). 

Fagures 2.5 through 2.8 show the calibration of the 
potentiometer settings including the frequency dependence 
for the coefficients used in the previous measurements. 
Figures 2.19 through 2.21 show the frequency response when 
the frequency compensated coefficients are used. 

The potentiometer settings were changed every 3 KHZ. 

During several of the "runs" it can be seen that there is 
MeieomeEshikt 1m the position of the experimental curve 
relative to the theoretical curve. This is due to 
uncontrollable drift of the basic clock frequency in the 
course of the datataking. 

The very latest set of data (Figure 2.21), was taken while 
adjusting a, and a, at every null, fieepuialbotenavine See said 


Ceefitcients according to the curves given in Figures 2.9 ~ 


Jans 


nn 





@eoneat that null. 
Table III shows the potentiometer-settings for a Ur Ss 


= = 1 a0 , ant One; iol = O67 15 b = Ora 7 





al 2 i 2 
aes ek 

3 KHZ 107 Z35 

6 " LOS .885 J 

o's 1.403 Nee 
i * ou .859 
iS Dy BOS 7 Tso .858 
18 =O" 862 





meer es 2.19 through 2.21 do not show any sign of the 
skewness of previous data sets, but the shape of the 
observed curves still deviate significantly from the 
theoretical curves, indicating that the proper coefficients 


have not been found for all frequencies. 
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FIGURE 2.9a. Theoretical and Observed Freauency Response of High Pass 


Filter; b, =-1.0, b, = 0.38 
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FIGURE 2,9b. Theoretic=1 and Ooserved Response of High Pass Filter; 
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FIGURE 2.10a. Theoretical and Observed response of Low Pass Filter; 


b, = 1,0, by = 0.38 
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FIGURe 2,10b. Theoretical and Observed Response of Low Pass Filter; 
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FIGURE 2,12b. Theoretical and Observed Response of Low Pass Filter; 
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FIGURE 2.15b. Theoretical and Observed Response of High Pass Filter; 
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FIGURE 2.17a. Theoretical and Observed Response of High Pass Piven. 
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FIGUEE 2.20a. theoretical and Observed Response of Low Pass Filter: 
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Wet eae ore LON TO MTIT—RADAR 


A. BOK Y 


ieee ctatement of the Radar Sienal Processing Requirement 


The majority of existing radars are pulsed radar 
systems. Thus the radar return is in effect a sampled signal, 
and the sampling width is the pulse width (PW) of the radar. 
Figure 3.la. is a time-amplitude plot of radar signals. 

At multiples of T= — the transmitted pulses are shown. 
The time interval between each transmitted pulse, called 
the Pulse Repetition Interval (PRI), target returns will 
appear. In Figure 3.la two stationary targets and one moving 
target are shown. The stationary target return may be either 
positive or negative depending upon the phase of the returned 
Signal. There is, however, no change from interval to interval, 
a fact which is exploited in the MTI canceller. 
The moving target however, is changing in amplitude every 
interval, at a rate which is the doppler frequency. 

Figure 3.1b shows how the pulse train of Figure 3.i1a 
would look as presented on an oscilloscope or A-scope. 
Here every PRI is shown on top of each other, and the moving 
target is exibiting the socalled "Butterfly effect" {6} 
While the stationary targets remain fixed in amplitude and 
time. 
To get the frequency representation of a periodic pulse 
[eati am thestime domain, cne can take the Fourier transform. 


The result would be a eumex |seype Weveshape. if only one 
Hs 
pulse were transmitted, the frequency spectrum would be 


Gontinuous, and the first null of the curve would appear at 
1 

PW’ But when transmitting a train of pulses spaced in time 
ee, = sa , ideally each target will return a square pulse 


Bee piecgqual to I. 


a 





The frequency spectrum of this periodic 


ion emean beetcuna inom the Fourier series if the pulse 
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where 
1 2 ~jnw,t 
Ce = 7 1 PGE) Myre dic »W4 = 
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CL js snown in Figure 3.2, where the sin x modulation 
o 
whose zero crossings with w- axis are determined by the 


pulsewidth. 


Pepe: For a radar with PRF of 800 HZ, a PW of 1.3usec 

and wavelenght » of 0.1 meter. 

The first zero crossing of the envelope will occur at 

769 KHZ. For a target radial speed ie. of 900 knots, 

using a formula from ref {6}; Doppler frequency; f,= 2v.,, 
yN 


d 





2 + 900 - 1852/3600 


5 7 = 9.26 KHZ 


i 


Thus it is obvious that the peaks of the repetitive 
"clutter spectrum" are essentially of equal amplitude in the 


frequency range of interest in MTI-radar (Figure 3.3). 


Pieuremsss 16 a small portion of Figure 3.2 (near the origin) 
and one moving target is introduced in this ideal frequency 


spectrum which consists of impulses. 
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Phau os 2. 


£ = 


Sin x 


envelope of many impulses at 


De lene 


Figure 3.4 shows the "real" clutterspectrum, i.e. the 


impulses have been spread out in frequency. There are 


several reasons for the nonideal shape shown: The energy 


in the returned echoes is spread out in frequence due to 


land target movement (vegetation swaying in the wind), sea 


waves moving, platform noise (own irregular movement), 


mG Poet SG On a finite RNUMber Of radar returns which stems 


from the fact that a radar antenna normally rotates. 


Thus we do not have an infinite sequence of target pulses. 


Rather, the number of 
iEypueaiiy, irom 10 to 
Figure 3.5 shows the 
antenna points at the 


Mmomeaqtiea the antenna 


"hits" on a particular target may be, 
50 per antenna rotation period. 

bursts of target echoes every time the 
target; the time between each burst 


rotation period. 


mmerewnes ¢.6 and 3./ the Fourier transform of a finite 


number of 6, respectively 10 hits are shown. The mathematics 


Ponenwteteda here. 


The fewer the number of hits, the more energy leaks into 
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FIGURE 3.4. A "Real" Clutter Spectrum 


frequencies where it may give raise to false alarm. 

Thus the finite number of hits not only leads to spreading 

of the clutter spectrum around multiples of PRF, but also 
reduces the maximum attainable cancellation ratio. For 

proper MTI operation, at least 10 hits per scan is needed. 
The Fourier transform of such a truncated series of target 
returns in which there is a moving target among heavy clutter 
return, is shown in Figure 3.8. Note that the component is 
also spread out from the ideal impulses, and the energy in 
the sidebands. 


b. Time Response Approach to MTI-Radar Analysis 
In this analysis we do not concern ourselves with 
the trailing edge of a typical echopulse, as our filter has 
no noticeable disadvantageous effect in this area. 
But since our filter is a recursive filter, the resirculation 
of signals lead to undesirable effects which must be 


investigated. 
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FIGURE 3.9. Delay Line Canceller 


Figure 3.9 shows the blockdiagram of a conventional MTI 
delay-line canceller {6}. Cancellation is obtained only when 
two succeeding pulses are equal in amplitude. The first and 
last pulse in a pulse train of a finite number of pulses will 
be passed by this filter. 

itmmerme recursive CCD-filter of Figure 2.3, the first 
pulse of such a fixed amplitude pulse train will be passed, 
the second pulse will not be fully cancelled. After the 
pulse train input is terminated, it will take a finite time, 
depending upon the values of the feedback coefficients ba 
and Dos before the filter output approaches zero. 
This time is measured in number of PRI's, and the effects in 
PPI presentation is that the target will spread out in 


Paice) the Girection Of the antenna rotation. 


Figures 3.10 and 3.j1 illustrates typical unit impulse 
responses of a first order recursive filter with feedback 
S@etticaents equal to 0.3 and 0.9, respectively. 

Example: A typical antenna rotation rate is 20 revolutions 


per minute, which means 3 seconds/scan or 8.33 msec/degree. 


For the case where by = 0.9 this means that the target 
would be extended 29/8.33 = 3.5 degrees beyond the actual 
position. This is considered to be a serious limitation, 
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FIGURE 3.10. Tame Response of First Order Recursive 


Filter, b, =O S 
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FIGURE 3.11. Time Response of First Order Recursive 


Filter, by = 0.9 
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emcaecOnsequentiy we have to keep this aspect in mind when 


designing filters of this type. 


Ds Applying uimemeem) Falter Hic eimenanie . 


The purpose of MTI -radar is to detect targets 
moving at or above a certain minimum speed, and suppress all 
echoes from stationary targets (clutter). 

To do this, one exploits the well known doppler effect, 


namely that the frequency of a returned echo from such a 


a 


target (moving), is- different from the transmitted frequency. 
Using detection” vechniques-in the recoiver {6} , {7},only 
“doppler frequencies are amplified and become the MTI- 
video which may be transmitted to a PPI for viewing. 

Thus we need a comb filter caestememe 2!) SS wWiion Can peject 


the clutter, and pass the doppler shifted signals. 











One disadvantage is that when the target's radial velocity 
corresponds to some multiple of the PRF, it will be 
Suppressed aiso (the target is then said to have a radial 
blind speed). To overcome this one must vary the PRI from 
Buese LO pulse, by PRF-jitter (randem variation of PRI), 

or PRF stagger (sequentially stepping the PRI between two 
values only). 

A proposed method for building a CCD + filter to facilitate 


such-asystem 1s outlined in chapter IV. 
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3. Bandwidth Requirements. 


The pulsewidth determines the bandwidth (BW) 


necessary to reproduce the signal. 


The employed delay lines have 96 stages each. The necessary 


clock frequency, fs is then: 


moore nRE = 800 HZ, f = lsowniaA, wihton Gc the effective 
sample frequency (f_). 


Using Nyquist's Sampling Theorem, fe = 2BW. 


The maximum bandwidth that can be faithfully reproduced is 


tmem 38.4 KHZ. This corresponds to PW = 13 usec, which is 


quite wide and unrealistic. Thus delay lines with considerable 


larger number of stages are needed. 


en EXPERIMENT 
ee Simulation of Radar Return Signal 


Since the size of the delay line was limited to 96 
stages it was decided to simulate the radar signal for which 
the delay line was able to respond. The cbjective was to 
provide radar return signals simulating fixed and moving 
targets with variable PRF and number of hits per scan. 

The simulation of PRF was implemented using 2 Beckman 9030 
function generator controlled by a sinewave generator 
(Wavetek Model 142). 

The function generator produces a square pulse whose lenght 
and amplitude can be varied for every period of the sinewave. 
Thus the PRF may be controlled by changing the Wavetek 
frequency (Figure 3.13). 
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Function Generator C 


~ FIGURE 3.13. Time Relationship Between Output from 


Funetion Generator and Sinewave Generator 


The pulse train thus obtained is the simulation of clutter, 


or fixed targets, since each pulse is equal in amplitude. 
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Figure 3.14 illustrates the simulation setup. A National 
Semiconductor MM 4616A switch (Figure 3.15) was used to 

perform the switching functions. The specifications of this 
arrangement is that each pulse through this switch had a 
different amplitude, simulating the moving target signal. 

On the oscilloscope, this signal should exhibit the "butter- 
fly effect" {6}. If needed, both stationary and moving target 
Signals could be fed into the filter simultaneously. 

The second switch was operated by another Wavetek M142, which 
was set such that the positive part of the squarewave was 
much shorter than the negative part. The positive part 

would turn the switch on, letting a finite number of pulses 
from switch one through. This was the simulation of the 
antenna sweep. However, the present simulation of this 
effect is limited because it has not been possible to 


Simulate the long time between two series of radar hits. 


2. Results. 


a. Stationary Targets with PRF = 800 Hz 





PecUiMmewstl/. Cancelling of “Stationary” Target, 


Mevizoental: 10 Wwsec/cem. Vert: .5V/cnm 


(2 





V =5V 
ee 






EPYromeeunction 


TO (ti) jer 


PECURE MoS. Anealogmowteen ochematic. 


TABLE IV 
Typical Radar Charachteristics 


Band: -band S-band {C-band X-band ‘-band 


| 
PPC QUCNG Ee) 2 ai | enoeniee— 12 5 PS uke 
(GHZ) 
Typical OO ele 500- 800- 3000- 
ee. 500 600 1000 | 3000 5000 







Pee ol bas 333- 200- 
ort MSEC 1250 usec 500 usec 
Maximum L62- 16.2- 
Theoretical 405 AO. 5 
Range n.m. 
Typical PW on 
0.5 usec 
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PrGuRE 3.18, Cancelling of Stationary Target. 
PREe=—eS00H4, 1oOrizemedl: 10 usec/cam. 





FIGURE 3.19. Cancelling of Narrow, Stationary Target. 
Pree mcUOHSertOorazeritalse1l0 USec/cm. 
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fcwneo.ce. Moving target Srenal-. PRF= 800 HZ 


Horizontal: 10 usec/em 


Figure 3.20 shows a "moving target" signal. The doppler freq- 
uency was at the center of first comb tooth of the filter. 
Pecan be seen that the combfilter did not cancel the signal. 


The insertion loss was quite small. 





picunt see, Moving and Stationary Targets. PRF= 800 HZ 


Horizontal: 10 usec/cm 
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CG. SOumeeisceTonarysamanMovineg Sipenals, PRF = 800 HZ 


Figure 3.21 shows an input consisting of both moving 
and stationary target signals. It can be seen that only the 


moving Signal is passed. 
d. Fixed and Moving Targets. PRF = 10 KHZ 


To try to eliminate some of the glitches, the effect- 
ive clock frequency was increased to 288 KHZ, so that the 
nulls of the comb filter appeared at multiples of 10 KHZ. 

The set of coefficients for the following results were 
(potentiometer settings): 


eecos85, a, = 0.975, a, = 0.87, b, = 0.5475, b,= 0.498. 


For all the pictures the vertical scale is 0.5 V/cm. 





PIGURBMOEe we ctattondry larcet meena. PRE = 10 KHZ. 


Hemazonte ) cca le molmuSec/ Cm. 


pugure 3.?2me@errespends to figure S218) >the difference 
being that now the PRF was 10 KHZ instead of 800 HZ. Note 
that the glitches were much smaller in amplitude at higher 


clock frequency. 


Tes 








PlGCUREMee coe OVang Target Ssipmal Only. PRF = 10 KHZ 


Horizontal scale: 50 usec/cem. 


e. Moving Target. PRF = 10 KHZ. 


Figure 3. 23 shows a moving target only, with 
doppler frequency, fa = 5 KHZ. Again it 1s noted that the 
moving target signal passed through the comb filter with low 


misertion loss. 





7, 
BiGuivemso.24. Movang Target, fi, Soe oie Ree 10 KHZ. 


Horizontal: 50 usec/em/ 
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f. Moving Target with Doppler Frequency Close to the 
Firsmewiee 
Figure 3.24 shows that the output amplitude indeed 
was attenuated when fa was close to the blind speed. 


In this case the first blind speed was the radial velocity 


of the target corresponding to fa =e 0! KHZ). 





PGURE 3225-0 MOvanp Target and stationary Target, 
fs =>) KHZ, PRE Seo Ke. 
Horizontal: 50 usec/cm. 


g. Stationary and Moving Targets, PRF = 10 KHZ 


Figure 3.25 shows that when the input consisted of 
both stationary and moving targets, only the moving target 


Signal was passed. In this case fs = 5 KHZ which was midband. 


h. Stationary and Moving Targets, Doppler Frequency at 
First Null 


Figure 3.26 depicts the same situation as in Figure 


3.25, only now f,= 10 KHZ which was the first blind speed. 


d 
Both targets were suppressed, and the glitches are seen. 
The output was almost identical to the output in the station- 


ary target case (Figure 3.22) 


ie 








FIGURE 3.26. Moving and Stationary Targets Superimposed. 


ta Set pars  Herpizontal: 50 
usec/cm. 


1. The Effect of a Finite Number of Stationary Target 
Signals. 





FIGURE 3.27. Truncated Series of Fixed Targets (15 hits) 


Horezontal.: OaSeamsec/cm. 


Beei@ees, 27 shows the input and output of the comb filter in 


the case where Switch 2 (Figure 3.14 ) was operated by a lew 
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frequency asymmetrical squarewave, to simulate finite number 
of "hits" on the target during one period of the antenna scan. 
The input to the switch was an infinite train of fixed 
amplitude pulses. The oscilloscope was synchronized to the 
Switch controlling waveform. Note that the first and last 
pulses in each truncated series of radar hits were passed by 
the comb filter. This was in agreement with the time analysis. 
The rest of the pulses were cancelled, except for the 
glitches, which appeared quite pronounced in this picture. 
Additional study is necessary to determine the cause of the 


glitches and their elimination. 





[ieURe a,2c. One Pulse from Figure 3.27. 


Homizontal scale: 10 ieec/cm. 


Figure 3.28 shows one of the individual pulses in the 
truncated trains of radar hits from Figure 3.27. The 
oscilloscope was now synchronized to the function generator 
(Figure 3.14). It is noted that the cancellation was poor, 
verifying the analysis of section III.A.1l.a, which predicted 
that energy would leak into sidebands when a finite number 


of hits were being processed, causing false alarm. 
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Mount o.29, Truncated Series of Fixed Targets (3 Hits). 


HeraZomeal: 20 sec en. 


Figure 3.29 shows a truncated series of only 3 hits. The 
oscilloscope was now synchronized to the squarewave generator 
controlling Switch 2, and the three hits cannot be seen 
Clearly on this picture for that reason. Note that the first 
two pulses were passed and the last pulse was also passed as 
well as several recirculated pulses, causing a smearing of 
the trailing edge. This phenomenon was desribed in section 
Witenes. >., and gs caused by recirculation of the signals. 
The Signal will die out gradually depending upon the size of 
the feedback coefficient. 


Figure 3.30 shows one of the individual pulses shown in 
Figure 3.28. Although this pulse was similar to Figure 3.28, 
it should be quite different. More pronounced false alarm 
Was expected in Figure 3.30 than in Figure 3.28. The reason 
for the poor result here was that the individual truncated 
series were too close, so that several consecutive series 
were superimposed on the scope. This was also the reason 

for the base line to appear even under the pulse in both 
Figures 3.28 and 3.30. 
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PGune 3.30. One pulse from Figure 3.29. 


Horizontal: bUMisec/cm. 


The relative brightness of the puise itself and the base-line 
under the pulse in Figure 3.28 indicated that there were more 
oscilloscope sweeps for pulses thanfor base line. 

But in Figure 3.30 the situation was reversed. This was 


reasonable, since here most of the time-axis was a base-line. 
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A. CONCLUSIONS OF EXPERIMENTS 


ieee EL iter Experiments 


faeeinaswbeen shown that the ClD-device can be used 


for implementing second order comb filters. 


ZeeeMiti=-Radar Simulation Experiment. 


Most of the simulation objectives were met except for 
miemsoimulation Of antenna rotation. Due to the limitation of 
instrumentation, the time between two series of radar hits 
is too short. A more realistic case would be to have a series 
of radar hits appear every two seconds. 

A major unresolved problem encountered in this study was the 
observed glitches. It must be solved before the sampled 
analog recursive CTD comb filter may be successfully employed 
in real radar systems. 

Characteristics of the present CTD-devices do not yet allow 
for enough bandwidth. More delay stages are needed. 

Until then only simulated results may be obtained, with 
rather restrictive ranges of pulsewidths and PRF's. 

However, deviations between theory and measured results have 
been noticed. It was found that both the delay line and its 
peripheral electronics do not have adequate bandwidths. 
Consequently, the coefficients used in experimental studies 
vary with frequency. It is believed that this is one of the 


major causes of deviation of experiment from theory. 


Pee COMMENDAT IONS 


eeteosen topres Vireetty Related to Present Thesis. 
In the filter study of chapter II, more work is 
needed to achieve implementation of frequency independent 
coefficients. Perhaps voltage controlled operational 


amplifiers should be used instead of potentiometers, and 
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also a frequency independent delay-line and its peripheral 
eumcuitry. 


In the MTI- simulation study, the problem is the glitches. 
To find the causes, and to eliminate them is important. 
Another area is to build a better simulation setup for 


examining the effect of finite number of returns. 


Peeeweposed Application to Staggered/Jittered PRF MTI- 
Radar. 


Not too many military MTI-radars today will be 
built that do not include either jitter PRF or: stagger PRF 
Capabilities. One of the major attractions of CTD-delay 
lines is the ease and speed with which the clock-frequency, 
and thereby the delay time, may be changed. Consequently, 
the null frequencies of the comb filter can be easily changed 
by changing the clock frequency. 
But there are major problems in applying the CTD-filter to 
such cases. This is because in jitter, or staggered PRF-radar, 
BHewERe Changes from pulse to pulse, and thus the clock 
frequency of the del) emma ue Demenamcedsrorrevery PRI. 
-It is seen that information which has already been fed into 
the delay--line (memory device), will be clocked out ata 
different clock frequency and then added in the output summer 
(Figure 2.3) with a signa employing a third PRI. The result 
would be meaningless. 
Figure 4.1 is a proposal to employ clocked CTD-delay lines 
in a jitter PRF tracking filter which could overcome this 
problem. It may solve the undesirable effects of long 
settling time associated with high feedback coefficients, 
eemogta-ed in section III.A.1.b. 
Pesnort explanation using difference equations, of this 


Careuit follows: 
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ist pulse arrives: 


Upper part of 
feelter effective: 


Video into lower 


em@eeurt, still 
tapping upper part: 


Sl opens. Video into 
mepe” part: 


ol closes, S2 opens, 


Video into lower part: 


Sl opens, S2 closes 


Video into upper part: 


¥(0) 


Y¥(1) 


eee (0) 


a X(1) + (bya ota, )xX(0) 


¥(2)sa X(2)+(b,a_ ta, )X(1) f 


1S = 


YC4)= 


yCS) 


% 2 


1 aX (0) +a 


4b ,X¢0) 


a Xx(3) + (bj, aot a, )x(2) 


aX(4) a (bj aot a,) KC 3) 


= aX (5) + (bia jta, )xC4) 


etc 


If the video train input is a train of fixed amplitude 


pulses, then: 


Y(n) + X(n){ (a (itb, ))ta, 3} 


There will be perfect cancellation if: 


Sacewie: Operation: 


a,@ + by? + a, = - 0 
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aa 1s a sawtooth waveform. A sample and hold circuit is used 
to establish a DC-yoltage level which determines the 
frequency output of the voltage controlled frequency genera- 
tors. Starting with frequency generator # 1, assuming that 
the video is being fed into the upper part of the filter, the 
binary counter will count 96 periods and then generate a 
trigger pulse, which controls the radar transmitter and all 
radar timing circuits and the various switches in the filter 
itself. The functions of this trigger pulse in the filter are 
to: 

- Reset the lower counter to zero. 

- Sample a to establish a new fA : 

=sowitch the video to the lower part of the filter. 

- Swatch off S2, the output of the lower delay line. 

~ Switch on Sl, the output of the upper delay line. 


Next, after another PRI and foo Counted tO $5, a tragger pulse 
appears (the counters will only produce a trigger pulse on the 
count of 96 after each reset). a is then set to a new value, 
counter # 1 is reset, the video is switched to the upper part 
mueethe filter, the switch Sl is opened and S2 is closed. 

Thus the direct signal will always be summed in one part of 
the filter with the delayed signal stored in the other part. 


This scheme may be extended to higher order filters. 
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‘ APPENDIX A 


FILTER DESIGN IN THE Z@ - DOMAIN 


It was desired to design a filter with 100 HZ across 
the null at an attenuation of at least 20 db down. 
From a pole-zero plot in the z-domain Figure A-1), the 
magnitude of the transfer function can be evaluated for every 
point on the unit circle, geometrically. 
By changing the pole and zero locations within the unit 
Seede sin a trial and error manner, and comparing the 
magnitude of the transfer function with the desired shape of 
the frequency response curve, it was possible to come fairly 


close after a rather short time: 





Results: Zeros: 245 =e Ocoee 7) Ont 


Poles: Poo = Ob, t a0.o 


This gives the transfer function: 


9.556 ~ 27 + 0.4556 277 


H(€z)= 
Ss woh ae =< 


Toe 2 0) ee 
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A computer simulation of’a filter using these cofficients 
is shown in Figure A-2. The desired width of the stopband, 
and attenuation was obtained. The filter was also built, 
and a few data points were plotted, but the implementation 


Mme cOetficilents were not correct. 
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